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We study theoretically the geometrical and temporal commensurability oscillations induced in 
the resistivity of 2D electrons in a perpendicular magnetic field by surface acoustic waves (SAWs). 
We show that there is a positive anisotropic dynamical classical contribution and an isotropic non- 
equilibrium quantum contribution to the resistivity. We describe how the commensurability oscilla- 
tions modulate the resonances in the SAW-induced resistivity at multiples of the cyclotron frequency. 
We study the effects of both short-range and long-range disorder on the resistivity corrections for 
both the classical and quantum non-equilibrium cases. We predict that the quantum correction 
will give rise to zero-resistance states with associated geometrical commensurability oscillations at 
large SAW amplitude for sufficiently large inelastic scattering times. These zero resistance states 
are qualitatively similar to those observed under microwave illumination, and their nature depends 
crucially on whether the disorder is short- or long-range. Finally, we discuss the implications of our 
results for current and future experiments on two dimensional electron gases. 

PACS numbers: 05.45.-a, 05.60.Cd, 72.50.+b,76.40.+b 
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I. INTRODUCTION 

Recent experiments on high mobility two-dimensional 
electron gases (2DEGs) have shown a variety of mag- 
netoresistance effects under intense microwave (MW) 
illuminationi2iMiSiSiLSi2ii£ One new phenomenon that 
has attracted considerable attention is MW-induced 
magnetoresistance oscillations. These are spaced accord- 
ing to the ratio of the MW frequency to the cyclotron 
frequency, showing a temporal commensurability between 
the MW field and the cyclotron motion. At high MW 
power the oscillations develop into zero resistance states 
(ZRS). 2 ' 3 The phenomenological theory of ZRS by An- 
dreev, Aleiner, and MillisAi has explained the effect to 
be a macroscopic manifestation of negative local resis- 
tivity imposed on the 2DEG by energy pumping from 
the microwave field. Several microscopic models have 
been propose d 12 ' 13 : 14 ' 15 to justify the periodic appear- 
ance of negative resistivity as a function of ratio of the 
MW frequency to the cyclotron frequency. The consen- 
sus is that the effect is of a quantum nature and related 
to the non-equilibrium occupation of Landau levels in a 
2DEG pumped by microwave radiation. 

A separate phenomenon involving geometrical com- 
mensurability oscillations has been extensively studied in 
2D semiconductor structures subjected to either surface 
acoustic wave (SAW) or static modulations ! 16 ' 17 ' 18 ' 1 ?' 21 ! 1 
Commensurability between the diameter of the cyclotron 
orbit, 2R C and the wavelength 27r/<7 of a coherent sound 
wave is known to manifest itself in magneto-oscillations 
of the dynamical conductivity a^ q of metals - the ge- 
ometrical resonance effect. This was used in studies 
of three dimensional metals to determine the shape of 
the Fermi surfaced Geometric resonance effects have 
also been observed in 2DEGs in the attenuation and 



renormalization of SAW velocity due to interactions with 
electrons 1 ^^ and in the drag effect due to SAWs. 22 
There has also been some theoretical study of these SAW- 
induced effects. 23 ' 24 ' 25 ' 26 ' 27 

In all of these studies, the parameter regimes used 
have been in the low frequency (u> — > 0) limit, such that 
the SAWs are essentially static. Experiments have re- 
cently begun to investigate the effect of SAWs on mag- 
netoresistance in a 2DEG4 6 " We have studied this effect 
theoretically, investigating both geometric and tempo- 
ral resonances. 28 For a spatially periodic SAW field, the 
commensurability effect can also be viewed as a resonant 
SAW interaction with collective excitations of 2D elec- 
trons at finite wavenumbersji 7 ^ 7 . enabling one to excite 
modes otherwise forbidden by Kohn's theorem^ These 
experimental developments require new theory for the re- 
sponse of a 2DEG to electric fields with both spatial and 
temporal modulation; this work was started in Ref. [2^ 
and is continued here. 

In this paper we study the non-linear dynamical effect 
in which SAWs induce changes in the magneto-resistivity 
of a high quality electron gas in the regime of classi- 
cally strong magnetic fields, lu c t 3> 1 and high temper- 
atures krjT 3> hu c (uj c — eB/m* is the cyclotron fre- 
quency, where B is the magnetic field, m* the electron 
effective mass, and r is the transport relaxation time). 
We have recently shown the existence of SAW-induced 
magnetoresistance oscillations 28 that reflect both tempo- 
ral and geometrical resonances in the SAW attenuation. 
There are two competing contributions to the resistivity 
corrections: one a classical SAW induced guiding centre 
drift of cyclotron orbits, the other a quantum contribu- 
tion arising from the modulation of the electron density 
of states. In Ref. |2^, we restricted the analysis to 2DEGs 
with short-range disorder (which leads to isotropic scat- 
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tering) , so as to understand the main qualitative features 
of the oscillations. We also restricted our analysis of the 
quantum corrections to oj < uj c . In the present work 
we extend our previous analysis of the quantum contri- 
bution to higher frequencies, u) > uj c , and address the 
experimentally relevant situation of long-range disorder 
which arises in modulation-doped systems and leads to 
small- angle scattering of electrons. 

The classical contribution to the resistivity correction 
originates from the SAW-induced guiding centre drift of 
the cyclotron orbits. For a SAW with frequency uj, and 
wavenumber q, propagating in the x direction with speed 
s = uj/q (we assume s is small compared to the Fermi 
velocity Vf), there is an anisotropic increase in the resis- 
tivity p xx (at high fields uj c t 1, this is equivalent to an 
increase of conductivity in the transverse direction a yy ), 
which oscillates as a function of inverse magnetic field. 
We show that at u> > uj c the resistivity change displays 
resonances at integer multiples of the cyclotron frequency 
uj « Nuj c . We find that the main difference between 
long- and short-range disorder is that long-range disor- 
der leads to an effective transport time r* = 2r/{qR c ) 2 
when qR c 3> lu/uj c . For small- angle scattering we ob- 
tain analytic formulae for the resistance correction for 
the limits uj c t* 3> 1 when (qR c ) 2 <C t/t s , and uj c t* <C 1 
for both (qR c ) 2 -C r/r s and (qR c ) 2 ^> t/t s , where r s is 
the total scattering time^ 

The quantum contribution arises from the modula- 
tion of the electron density of states (DOS), 7(e), and 
consequently, from the energy dependence of the non- 
equilibrium population of excited electron states caused 
by Landau level quantization. We follow the idea pro- 
posed in Ref. [lj to explain the formation of ZRS^iil 
under microwave irradiation with oj > u) c . We show 
that in the frequency range r" 1 < uj < uj c the quan- 
tum contribution suppresses resistivity both in p xx and 
p yy and persists up to temperatures fc^T ^> Hoj c and fill- 
ing factors v 1 where no Shubnikov-de Haas (SdH) 
oscillations would be seen in the linear-response conduc- 
tivity. We propose a new class of ZRS, in which ge- 
ometric commensurability oscillations overlay the ZRS 
that would be found in the microwave (qR c — > 0) limit 
for a short-range potential. For a long-range potential 
there are ZRS linked to geometric commensurability os- 
cillations, which are enhanced (by 0([qR c ] 2 )) over those 
induced by isotropic scattering. 

The paper is organised as follows. In Sec. [n] we give 
qualitative arguments to determine the form of the clas- 
sical magnetoresistance oscillations in the presence of 
short- and long-range disorder potentials. In Sec. IIIII 
we obtain these results rigorously using the classical ki- 
netic equation, and discuss screening of the SAW field. In 
Scc. lIVI we give our analysis of the quantum kinetic equa- 
tion at both low and high frequencies. Finally, in Sec. 
we discuss our results, in particular their implications for 
experiments. 



II. QUALITATIVE ANALYSIS OF CLASSICAL 
MAGNETORESISTANCE 

It was shown by Beenakker— that the magnetoresis- 
tance oscillations in a 2DEG in the presence of a spa- 
tially modulated electric field can be understood from a 
semi-classical point of view by considering the guiding 
centre (E x B) drift of cyclotron orbits which lead to 
enhanced diffusion. We apply a similar method below 
to calculate the guiding center drift in the presence of 
an electric field with both temporal and spatial oscilla- 
tions, for both isotropic scattering (short-range disorder) 
and small- angle scattering (long-range disorder). Note 
that the following qualitative analysis gives quantitative 
results applicable only to the high-field regime, uj c t I. 

The dynamical classical resistivity change 8 c p xx can be 
tracked back to the SAW induced drift Y(t) (along the y- 
axis) of the guiding center (X, Y) of an electron cyclotron 
orbiliiS and the resulting enhancement of the transverse 
(y-) component of the electron diffusion coefficient, D yy , 

S c p X x _ 5 D yy 
Pa D 

where p = 2/je 2 VpT is the Drude resistivity and D n = 
R 2 /2t is the unperturbed diffusion coefficient (where the 
cyclotron radius R c = vp/a> c ). The drift is caused by an 
electric field E^q cos(qx — wi)x, where x(t) is the position 
of the particle. To lowest order in E uq , the contribution 
of the SAW to the guiding centre drift velocity is 

Y(t) ~ ^CQB[qx(t)-ut], (I) 

X(t) = 0, (2) 

where x(t) is the position of the particle neglecting the 
effects of the SAW, but including the effects of disorder 
scattering. The change in the electron diffusion coeffi- 
cient due to the SAW is then 

SD yy = £ (Y(t)Y(O)) dt, (3) 

where we average over all particle trajectories x(t) in the 
disorder potential. 

A. Short-range disorder: isotropic scattering 

For isotropic scattering, the particle performs free cy- 
clotron orbits, x(t) — R c sm(uj c t + <po) + Xq, up until 
a scattering event, after which the subsequent motion 
in the SAW potential is uncorrelated with its preceding 
motion, provided qR c ^ I . In this case, averaging over 
the trajectories and scattering events, the change in the 
diffusion coefficient due to the SAW is 

6D m =f dte^ fj J* !£d* nt)m (4) 
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where ip = qXo, and Y(t) is calculated from Eq. us- 
ing the free cyclotron motion. Using the Bessel function 
identity 



E 

N=-o 



J N {z) 



JN<f> 



(5) 



one can obtain the frequency and wave number depen- 
dence of the magnetoresistance effect: 



S c p x 



5D V 



Po 



where 



Dn 



N=-o 



J 2 N (qR c ) 



l + (w- 



SAW 



f» 77?, 



(6) 



(7) 



is the dimensionless SAW amplitude. The electric field 
E^q discussed in this section is a travelling wave, rather 
than the standing wave situation discussed in Sees. II I II 
IIVI and Ivl It is related to the E uq discussed in Sec. IIIII 
by E uq — 2E ulq , implying £ = 2£. If we include another 
travelling wave to generate a standing wave, then the re- 
sult in Eq. © should be multiplied by a factor of 2, and 
this reproduces Eq. ®. In Eqs. © and 0, E^ w is 
the SAW longitudinal electric field, a scr = x/27re 2 7 the 
2D screening radius, tp the Fermi energy, I — vpr the 
mean free path, x the background dielectric constant, 
and 7 = m/irh 2 the electron density of states. Equa- 
tion |JBJ includes the Thomas- Fermi screening of the SAW 
field by 2D electrons, E uq — qa scl -E^ w , which we dis- 
cuss in detail in Sec. lIII ATI At large qR c we can expand 
the Bessel functions to get 



5 c Px 



■nqRc 



E 

JV=-oo 



cos 2 (qRc + jgr-_f) 
1 + (lu- Nlo c ) 2 t 2 



(8) 

From Eqs. JJJJ and (JBJ we can see that there are a se- 
quence of resonances at integer multiples of the cyclotron 
frequency, lo « Nu> c . The widths of these resonances are 
controlled by lot, with large values of lot leading to very 
narrow resonances. The oscillations for even harmonics 
are in phase with the Weiss oscillations of the static po- 
tential, whilst those of odd harmonics are tt out of phase. 
This can be understood by noting that the main contri- 
butions to the drift occur when the electrons are moving 
parallel to equipotential lines. For odd harmonics the 
phase of the potential at the half-orbit point is opposite 
to that for a static potential, and hence the cancellation 
and reinforcement effects that lead to minima and max- 
ima in the resistance are interchanged between the static 
and dynamic cases. This is illustrated via a comparison 
of the two situations in Fig. ^ Alternatively, this can 
be seen from examination of Eq. © , since the static case 
(lo = 0) is dominated by the N = term in the sum, 
whereas all N contribute in the dynamic case. 



Even Harmonic 



t = 



t _ 7t | Force due to static potential 

to 

| ^ T Force due to dynamic potential 



Odd Harmonic 



t 



FIG. 1: Comparison of electron motion in a static and a dy- 
namic periodic potential. When the forces on opposite sides 
of the orbit are in the same direction, there is a maximum 
in the magnetoresistance. The horizontal lines correspond to 
maximum positive (solid lines) and maximum negative (dot- 
dashed lines) values of the field. In the dynamic case the 
direction of the field oscillates with frequency uj. Thus, when 
u) is an odd multiple of the cyclotron frequency lo c , there is 
an interchange of the resistance maxima and minima as com- 
pared to a static potential. 



In the regime lo c t 3> vjr/s, the resonances are very nar- 
row and appear to display a random sequence of heights, 
rather than the linear in B dependence evident in Figs. [21 
and |31 reflecting the dependence on the geometric reso- 
nance conditions. In the intermediate frequency domain, 
a natural regime for GaAs structures 



T- 1 < LO < LO. 



with densities n e > 10 cm at sufficiently high mag- 
netic fields, the classical oscillations take the form 



S c Pxx 
Po 



4.s 



^£ 2 4(qR c ). 



(9) 



The competition between electron screening effects (oc 
q 2 ), the dynamical suppression of commensurability by 
the SAW motion (oc lo~ 2 ), and the relation uj/q = s, 
means that the form of these oscillations is independent 
of the absolute value of the SAW frequency, provided 
that conditions r _1 < lo and Vf > s are satisfied. Note 
however, that for a fixed SAW amplitude, there is still 
q dependence of the SAW field E^ w since the electric 
field induced in the 2DEG by the piezoelectric coupling 
is a function of qd, where d is the distance between the 
surface and the 2DEG^i 

We illustrate the interplay between dynamical res- 
onances in the time and space domains by plotting 
Eq. © as a function of lo c /lo for the following experi- 
mentally relevant parameters. In GaAs based 2DEGs, 
to* = 0.068 TO e , (where m e is the electron mass) and 
s ~ 3000 ms _1 . The highest SAW frequencies that have 
been used in experiments on 2DEGs are ~ 10 GHz^ we 
consider this fequency w ith t he sample densities and mo- 
bilities reported in Refs. 1 3. In Ref. 01 the mobility fi is 
3 x 10 6 cm 2 /Vs, and the density n e — 2 x 10 n cm~ 2 , 
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FIG. 2: Magnetoresistance for a short-range potential using 
the parameters of Ref. 1, and assuming lo — 2n x 10 GHz 
[lot = 7) and vf/s = 63. We choose a dimensionless SAW 
amplitude of £ = 0.01. 



corresponding to lot = 7 and vjr/s = 63 (note that 
vp/s = qR c when u> — lo c ). The magnetoresistance 
for these parameters is plotted in Fig. [21 In Ref. 0, 
(j, = 2.5 x 10 7 cm 2 /Vs and n e = 3.5 x 10 11 cm -2 , making 
it one of the highest mobility samples yet fabricated. For 
these parameter values, lot = 60, and vf/s = 84, and 
we find a magnetoresistance trace as shown in Fig. [3] At 
lower frequencies in such a high quality 2DEG, we ex- 
pect the width of the resonances to broaden and become 
similar to those shown in Fig. [21 




FIG. 3: Magnetoresistance for a short-range potential using 
the parameters of Ref. 2, and assuming lo — 2tt x 10 GHz, 
which implies lot — 60 and vf/s — 84. We choose £ = 0.001, 
a factor of 10 less than in Fig. 2. 

The dynamical mechanism just described dominates 
in a classical electron gas. No redistribution of electron 
kinetic energy (due to SAW absorption) will additionally 
change the magnetoresistance until the 2DEG is heated 
to a temperature kgT e > ftw c kp/q where geometrical 
oscillations become smeared. The essential assumption 



leading to this statement is that electron single-particle 
parameters (velocity, v, and t -1 ) vary slowly with energy 
at the scales comparable to the Fermi energy and can thus 
be approximated by constants. However, for high-quality 
2DEGs (as formed in modulation doped GaAs devices) 
the typical disorder potential is not simply short-range 
(as assumed for the isotropic scattering model considered 
above), but is dominated by a long-range part due to 
the coulombic potentials associated with remote dopants; 
this long-range disorder leads to small-angle scattering. 
We now turn to consider this situation. 



B. Long-range random potential — small angle 
scattering 

To calculate the resistivity correction in the presence of 
long range disorder, we consider an electron undergoing 
cyclotron motion and subject to random (small-angle) 
changes in direction. We express the position of the par- 
ticle as 



r(t) = R(t) + -Rc(sin </>(*), cos <j>(t)), 



(10) 



where R(i) is the position of the guiding centre and <j>(t) 
is the angular position around the cyclotron orbit. In the 
absence of disorder, <fi(t) = 4>q + io c t and R(t) is constant 
in time. The effect of a disorder potential on <j>(t) may 
be represented by random changes: the scattering events 
are separated by a characteristic time r s (sometimes re- 
ferred to as the quantum lifetime), and at each scattering 
event <f> jumps through a (random) angle of magnitude 
5<j), which is related to the lengthscale of the disorder 
potential L by 6(f) ~ l/(kpL) where kp is the Fermi mo- 
mentum. At any scattering event only the direction of 
motion changes; the position r(i) is constant, so Eq. (|1(J|> 
implies that for an instantaneous change 5<j> in </>, there 
is a change in guiding centre position of 



SR = R c (— cos (j), sin (f>)Scj). 



(11) 



In this section we consider the limit of a long range po- 
tential, for which 5(f) <C 1. In this case the scattering of <p 
can be viewed as leading to a continuous diffusive motion 
with a diffusion constant (in angle) of order ~ (<5</>) 2 /r s . 
Similarly, we can view the guiding centre co-ordinates 
as undergoing continuous diffusive motion, provided the 
typical jump [Eq. in guiding centre co-ordinate is 

small compared to other relevant lengthscales, in partic- 
ular the wavelength of the SAW, R c 5(j) < l/q. The dif- 
fusion constant for the guiding centre is ~ i? 2 (5</>) 2 /r s . 
Note that the diffusion constants for <f> and R depend 
only on the transport relaxation time r ~ r s /(<5</>) 2 , so 
only r enters the theory. In terms of r, the conditions 
for validity of this theory (S(f> <C 1, R C S4> < 1/?) can be 
written t/t s ^> 1 and (qR c ) 2 <C t/t s . 

To model the continuous diffusion of <f> and R, we in- 
troduce two sources of Gaussian noise, C\ and £2, with 
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we find a resistance correction of 



(£ i (t')£ i (t)} = T i 6(t-t'), 



(12) 



(where i = 1,2) and write the effects of scattering on the 
phase of the orbit and on the ^-component of the guiding 
center position as 



4>{t) = 4> + u c t+ / Ci(t')dt', 
Jo 



(13) 

X(t) = X + R c I £ 2 (t')dt'. (14) 
Jo 

The two sources of noise are related through Eq. (|llf> 
Ci{t) = -R c Ci(t) cos cj>(t), 

which requires that Ti = 2T 2 in the situation of interest 
where we average over all trajectories (hence over </>o). 
From a direct evaluation of the classical Kubo formula 
for the diffusion constant for this model in the absence 
of a SAW potential, we find that 



Ti = 2r 2 = 



(15) 



where r is the conventional definition of the transport 
relaxation time, in terms of which the Drude expression 
for diffusion constant is D xx = D yy = hvpr/(l + co 2 t 2 ). 

We simplify subsequent calculations by ignoring cor- 
relations of Ci and C 2 beyond those in Eq. I|15l) . which 
we expect to be accurate for lu c t* 3> 1, (with r* defined 
in Eq. (|20|l below) in which case the particle is able to 
explore all values of <fi(t) on the timescale of the rele- 
vant scattering time. We can then calculate the change 
in the diffusion coefficient from the guiding centre drift 
as in Sec. Ill Al for isotropic scattering, averaging over the 
electron trajectories 



SD yy = 



dt((Y(t)Y(p) 



F 2 
B 2 



dt 



Ci,C 2 

d<j)o 
2tt 



jiqR c (sin <f>(t)— sin 0o) 



x e 



iqR c J* C2(t')dt / — iujt 



Cx,Ci 



(16) 



where our assumption that C\ and C 2 are uncorrelated 
allows us to perform the average over noise: 



x ( e 



SDyy oc E 2 q Re 



Using 



^iqRc (sin <f>(t)— sin t^o) 



JqRcf* C 2 (t')dt' 



1 ™- ' ! dte~ wt 

Jo 

iqR. c (sin 4>{t) — sin 



2tt 

iqR c p C 2 (t')dt' 



(17) 



Ci 



= ^ JN(qR c )JM(qRc) 



N,M=-c 



x e 



i(N-M)4> 0f ,iNu c t p -2-^±t 



-r 2 t 



Po 



1 2 

4 S 2 t 



E 



JN(qRc) 2 T Nq 



N=-c 



1 + {NUJ C - LU) 



'Nq 



where 



TNq 



N 2 



(qRc 



(18) 



(19) 



The form of this correction is that found for isotropic 
scattering, but with rjv g replacing r in the correction to 
the diffusion coefficient,^ The result in Eq. 1)18(1 can be 
related to the result found in Eq. (|58|l in a similar manner 
to the case of short-range scattering. 

In the experimentally relevant range of parameters, 
lo/lo c <C qR c , and noting that in the vicinity of reso- 
nance N, uj/uJc — N, we define a scattering time r* <C t 
as 



„2 c>2 

q R c 

It 



(20) 



and tn q ~ t* in our regime of interest. This quantity 
enters in our discussion of the kinetic equation for small- 
angle scattering in Sec. IIII Bl Note that the use of this 
result is also subject to the constraint uj c t 3> 1 with the 
additional requirement qR 2 <C I, which is discussed in 
more detail in Sec. IIIII 



III. CLASSICAL KINETIC EQUATION 

We analyze the classical kinetic equation for a 2DEG at 
temperature fcsT < hu} c kp/q irradiated by SAWs. Our 
approach is to solve the kinetic equation for the electron 
distribution function 



f(t, x, ^e) = f T + y] e — J2 



(21) 



using the method of successive approximations. Here, 
/t(c) is the homogeneous equilibrium Fermi function, 
and the angle ip and kinetic energy e parametrize the elec- 
tron state in momentum space. Each component f m de- 
scribes the rn-th angular harmonic of the time- and space- 
dependent non-equilibrium distribution. To describe lo- 
cal values of the electron current and the accumulated 
charge density, we use the energy-integrated functions, 
9uq = Io°° dcf™ q - The relaxation of the local non- 
equilibrium distribution towards a Fermi function char- 
acterized by the value of local Fermi energy, €p(t, x) (de- 
termined by the local electron density n(t, x) oc g (t, x) — 

f™def°(t,x), where f°(t,x) = J ^f(t,x)) and the ki- 
netic equation is 



£f = Cf, 



(22) 
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where 

C = dt + v cos ip d x 



eE . 

uj c sin p 

P 



d v + evE cos <pd e , 
(23) 



with C the collision integral, E the electric field and p 
the electron momentum. 



A. Isotropic scattering 

In the presence of a short-range potential, the scat- 
tering is isotropic and we can use the relaxation time 
approximation for the collision integral, which is 



Cf 



f-f° f°-f T (e-e F (t,x)) 



(24) 



where we distinguish between the elastic scattering rate 
t -1 and energy relaxation rate . 

The dynamical perturbation of the distribution func- 
tion can be found from time/space Fourier harmonics of 
Eq. 11221) at the frequency /wave number of the SAW, 



1 w 

cL H % V iqR c cos <p 

UJ C T LU C 



f m = *ft>), (25) 



K 



2tt 



dip f v dip 



2tt 



G{ip,0) 



= E T 



J 2 N (qR c ) 



+ ir(Nuj c — u>) 



(28) 



which has the properties 

K(w, q) = K*(-w, q) = K(w, -q), 
and is obtained from Eqs. Il2(ill and H27|) using Eq. (JSJ. 

1. Screening and Dispersive Resonance Shift 

The electric field E in Eq. l)2.'{[l is the combination of a 
homogeneous DC field Eoo and the screened electric field 
of the SAW, E(t, x) = ^ uq E^e^'^x, found from the 
unscreened SAW field via E uq = E^ w /k(uj, q), where 
K(u),q) is the dielectric function of the whole 2D struc- 
ture. The density modulation n wg = 7<?^ ? induced by 
the SAW sets up an induced field E^ = -^f^TsSg, 
that we take into account at the level of Thomas-Fermi 
screening, so that E uq = E^ w + E™^. In the analysis of 
screening, the DC part of the electric field can be ignored 
(/oo = 0), and self-consistency yields 



*(¥>) = 



(deft 



v cos ip d e 



simp 
P 



-f 



Woo + /r), 



where we include the unknown perturbation of the 
time/space averaged function, /oo, related to the DC cur- 
rent to lowest order in E uq . We note that as Ti n only con- 
tributes to heating of the 2DEG, we can ignore it here for 
the AC part of the distribution function, but must retain 
it for analysis of the DC part of the distribution func- 
tion. In our expression for \E' we neglect contributions 
from higher Fourier harmonics of /, such as ftu>2q, /2uO, 
and /o2 9 : since these only affect DC transport at quartic 
order in the SAW field, whereas our interest is only in 
effects that are quadratic in the SAW field (i.e. linear in 
the SAW power). We also omit terms involving E 00 ft q , 
since they contribute to DC transport only beyond the 
linear response regime in the DC field. Equation 1)25(1 can 
be formally solved using the Green function G(ip, ip) 



9u)q 



eE^q 1 — (1 — iujr)K 
~iq l-K ' 



and 



«(^,t?) = 1 



1 l-(l-iojT)K 



9 



1 -K 



(29) 



(30) 



In the limit that qR c 3> 1, which is our regime of interest, 
the dielectric function becomes 



K(u),q) = 1 + 



1 - (1 - iurr) K 



(31) 



and for most SAWs, a scr |g| <C 1, which implies k ~ 
l/(GscrM)- When we account for dispersion, we find that 
the system has resonances at w = Nlo c + A jy where A jy 
is the dispersive shift of the N th resonance. 

We use Eq. (|3(J[1 to find the eigenmodes of the system. 
Setting k = and inserting the expression for K, one 
finds to leading order in lu c t, 



fu>q{<P) = 



G(ip,<p)V(<p)dip, 



G(cp, tp) 



(26) 
(27) 



which allows for an infinite range of variation of ip whilst 
guaranteeing periodicity of the solution f uq (<p). We also 
introduce the useful quantity 



where 



A = NLO c J 2 N {qR c ) 
N a scl .\q\ + l-J 2 N (qR c )~NA N : 



a n = 2^ 

p=l 



p 



(32) 



(33) 
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In the limit qR c — > we recover the magneto-plasmon 
dispersion, Ax ~ \q\R^u) c /a SCI , whereas in the limit of 
qR c ^> 1 screening becomes important and Ajv -C Nu> c . 
The crossover between the two regimes, which may be 
evaluated by considering the denominator of Eq. (|32ll . 
occurs at a wavenumber qoR c ~ a scr /R c <C 1. In our 
discussion below we take into account the screening of 
the SAW field. 



2. Magnetoresistance oscillations 

To find the steady state current, we analyze the 
time/space average of the kinetic equation in Eq. (2D) 
and take into account the dynamical perturbation ft q 



T- 1 - T- 1 



loo 



in j*0 



ev • E, 



on 



« sm( P* 

V COS If O e Oo; 

P 



deft 



ft q - (34) 



In our analysis of classical magnetoresistance oscilla- 
tions we assume that Ti n 3> t and since we are not 
interested in the heating associated with this term we 
drop it here. However it is important for our analy- 
sis of the first quantum correction to the classical re- 
sult that we present in Sec. IIVI We substitute the so- 
lution Eq. into Eq. (J33}, keeping track of the ef- 
fect of the perturbation of the time/space averaged func- 
tion /go on f uq . This procedure automatically includes 
SAW-induced non-linear effects. We multiply Eq. 
by (2aj c /ev F )ve~ lip , integrate with respect to e and cp, 
then use the relation between the x and y components of 
the DC current, j x — ij y = e^vpg^Q and the harmonic g^o 
(note that electrical neutrality requires g^ = 0)> which 



gives 



2u n 



ev. 



2tt 



UJ C T 



v cos ip d e 



v/oo 
smtp 



P 



fwq 



Jx 



v|re 2 7/2 2 



In 



leE, 



K 



1 -K 



iiv c T + 1 

+ I — — ~ I [Jx ~ IJyl 



v|re 2 7 /2 



IPX 



^00- 



(35) 



is the Drude resistivity tensor. Thus we find the resistiv- 
ity corrections 



o~ c Px 



Po 

S C Pyy 

S C P V X 



ujq 



elE, 



CF 



Re 



K 



1 - K 



(36) 



Po 



5 C Pxy 
PO 



2 ^ 

ujq 



leE u _ 



(F 



Im 



K 



1 -K 



and with the use of E m = E^ w / K (uj,q) ~ qa scr E^ w , 
Eq. (|28p. uj c t 3> 1 and qR c > 1 we formally justify the 
result in Eq. (jSJ. The magnetoresistance correction is: 



° Pxx = 2 (qlSf ^ 



Po 



^ 1 + (w - Nuj c y r 2 



(37) 



Strong damping (uj c t < lj 

At low magnetic fields, there is experimental 35 and 
theoretical 3 ^ 3 ^ evidence that there is exponential damp- 
ing of Weiss oscillations, and it seems natural that similar 
behaviour should be observed for SAW-induced oscilla- 
tions. We explore this question and find the functional 
form of the damping for isotropic scattering in this sec- 
tion, and for small angle scattering in Sec. IIII Bl In the 
strong damping limit (lu c t < 1), we investigate Eq. I|28|l 
when qR c 3> 1 and find the values of ip and <p that lead 
to a stationary phase. We then integrate over fluctua- 
tions about each point of stationary phase. These saddle 
points are ip — tt/2 and (p — 37r/2, and <p takes values 
which are any positive or negative odd integer multiply- 
ing 7r/2 such that <p < (p. When we sum the results of 
integrating about each saddle point, we get (to lowest 
order in e~ "<= T ) 



K = 



1 

LU C T 



1 



\qR c \ qR c 



sin(2qR c ) 



(38) 



In the limit that the magnetic field goes to zero, with 
ql ^> 1, this leads to a resistance change 



S c Px 



Po 



2ql£ 2 , 



(39) 



B=0 



and if we consider the resistance change 5p xx , we find it 
is 



and this can be used to determine the classical SAW in- 
duced change of the resistivity tensor, S c p. The relation 
between the electric field and current is E = pj + <5 c pj, 
where 



Po 



1 LU C T 
-UJ r T 1 



k 1 o2 



S C p X x = 

7re z pf 



e 2 



(40) 



which is independent of disorder. Summation over ±u> 
and taking the imaginary part of K as in Eq. (|36|l leads 
to the same result that 5 c p xy = 0. 
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We interpret this as a SAW-induced backscattering 
contribution to the resistance, which will dominate in 
the limit that the SAW wavelength is much less than the 
mean free path. We ignore the contribution to the resis- 
tance from channeled orbits, which can also lead to a pos- 
itive contribution to the magnetoresistance in the small 
magnetic field limit. 37 - 3 & The condition for the existence 
of these orbits is that the force from the screened SAW 
field is larger than the Lorentz force, i.e. £ > ^S-M^ 
and we assume that £ is sufficiently small for their contri- 
bution to be ignored (this is generally the case over most 
of the magnetic field range that we show in our figures). 



B. Small angle scattering 

A long range disorder potential leads to a non-isotropic 
scattering probability, and this implies that there are two 
scattering times that we need to take into account. One 
is the total scattering time, t s ^ and the other is the 
momentum scattering time, r, and in GaAs heterostruc- 
tures, r ^ r s . In the limit that (qR c ) 2 <C — , the disorder 
potential leads to diffusion in angle, and can be studied 
by replacing the collision integral in the kinetic equation 
by a term involving two if derivatives^ The dynamical 
perturbation of the distribution function, Eq. 125(1 is thus 
modified to read 



d v - i — 

L0 r 



iqR c cos if 



where 



eE_ 



v cos if d e 



smip 



0_: 



f ug = (41) 



(/oo + /t). (42) 



In Eq. ((41(1 . small angle scattering is introduced in the 
form of diffusion along the Fermi surface and is taken 
into account by the term ■ Now, if we let f uq = 

h wq {ip)e' iqR " s{v ' v , and substitute this into Eq. (gTJ), then 
solve in the limits that qR c /oj c T <C 1 and u}/uj c *C qR c 
(the second condition allows us to ignore the term d^h), 
we can solve the kinetic equation as before to get 



d<fG(ip,<f)V(tf), 



(43) 



where 



G(f,f) 



exp iqR c (sin dp — sin if) 

1 lo\ 
~ i— )(& ~ <P) 



1 



2w c r 4 



(sn^-sh^H , (44) 



and r* was defined in Eq. ((20(1 . Our results for isotropic 
scattering are modified by replacing G by G. The long- 
range potential problem is then reduced to a calculation 



of K for the modified Green function, which we define in 
analogy with Eq. 128|) as 



(i 



2?r 



(45) 



We can write the following exact expression for K: 
J s+2m-2p{.qRc)J s{qR c ) 



K 



s,m,p— — oo 



+ i((s + 2m)uj c — ui)t* 



' ^ In f^-j I , (46) 



2uj c t 



2u r T* 



where I m (x) is the modified Bessel function of the first 
kind. We study K in the weak damping (lo c t* 3> 1) and 
the strong damping (lu c t* <C 1) limits. 



Weak damping (lu c t* ^> 1) 

In the weak damping limit, we can make use of the 
asymptotic expansion of the modified Bessel functions 
for small argument, and need only retain the m = p = 
terms. K has the same form as K, except that r is 
replaced by r*, i.e. 



K = 



N 



(qRc) 



JV=-o 



i(Nu c 



(47) 



where \ 



was introduced in Eq. I(20() . 



Strong damping (ui c t* < lj 

In the strong damping limit when qR c 3> 1 and 
(qR c ) 2 «C -p, we investigate Eq. (|4*51) and use a simi- 
lar saddle point procedure to the one we used for strong 
damping in the case of isotropic scattering. After adding 
the contributions from integrating about each saddle 
point, we get (to lowest order in e~ ^>= T ' ) 



K 



UJcT" 



\qR c \ qR c 



sm(2qR c ) 



(48) 



Note that this is the same form as the expression for 
K for isotropic scattering in the limit w c r <ti 1, except 
that r replaces r* . If we want to investigate the limit in 
which the magnetic field goes to zero, then the kinetic 
equation [Eq. I(41(l ] as written previously is inapplicable 
when (qR c ) 2 3> 



30 



In our calculation above, there is 



F , which arises naturally 
as a result of our saddle point analysis. In Ref. 30 an 



a damping factor of D = e" 
as a result of our saddle pi 
alternative approach was used to calculate the damping 
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factor in the low magnetic field limit. In that approach 
the damping factor above is the high field limit of 



D\ = exp ■ 



(49) 



where r s is the total relaxation rate. If, as in Sec. Ill Bl we 
assume that phase and guiding center corrections are un- 
corrected, we can calculate a second damping factor as- 
sociated with phase corrections in addition to the guiding 
center contribution in Eq. I|49|) . using the same method. 
This gives 



which is very similar to Eq. i|31|) with K replacing K. 
When qR c 3> 1, screening of SAWs is identical for both 
short and long range potentials. Analysis of Ajv and An 
as in Sec. IIII ATI leads to the same qo at which screening 
becomes important. 



2. Magnetoresistance oscillations 

In this section we derive the magnetoresistance oscil- 
lations analogously to Sec. IIII A 21 The system of equa- 
tions that we wish to solve is 



D 2 = exp 



JCTS l + (^)\ 



(50) 



and at moderate fields, we may approximate D = D\Di, 
so that 



K 



LO c T" 



1 2 *™ , 

—— + —De »c sm(2qR c 
\qRc\ qR c 



(51) 



If we assume that the two effects are correlated (as 
we should in the (qR c ) 2 S> — limit), then we get the 
expression 



D = exp 



1 



(ql sin0 — lot) 



u c t s Jo {qlsiiKp - lot) 2 + (lo c t)(lo c t s ) 

which implies that in the lo c t — ► limit, with ql 3> lot 3> 
1 and — > 1 



D ~ exp 



1 



UJT I T s 

-^UOrT, 



(53) 



which gives lim^^o D = e u '" T = ^ lim CJt ,_,o ~D\Di — 
1. Screening 

In the above discussion, the solutions obtained for f uq 
in Sec. IIII Al and for f uq when there is small-angle scatter- 
ing differ in that the kinetic equation for isotropic scat- 
tering contained the term g® q /uj c T, which is absent here. 
When we solve for g® q , and hence the dielectric function, 
Eq. (|3L)[) . in the presence of small-angle scattering, we get 



fJ. 



iq 



1 - (1 - iLOT*)K 



which implies a dielectric function 

1-{1-IUT*)K 

k{u, q) = H n , 



(54) 



(55) 



V COS ipO e Ou. 



-OeJT H C Jo, 



L0 C 



(56) 



in combination with Eqs. (|41|l and (|42|l . The solution 
of Eq. |gTJ for is shown in Eq. g3J). We take the 
equation for / , and integrate with respect to energy and 
tp after multiplying by ve~ l<p , as before, and obtain 



6E. 



oo 



7 A R y 



2w r 



em 2 v F 



eE, 



to, 
dtp 
2^ 



dtpG(ip, 0), 



which is the same as we found for isotropic scattering, 
except that G replaces G. The resistivity correction is 
thus 



Po 



elE t , 



CF 



-Re 



{*} 



(57) 



As in the case of short-range scattering there are no cor- 
rections to any other components of the resistivity tensor. 
Unlike the case of isotropic scattering there is not a fac- 
tor of 1 — K in the denominator of Eq. I|57() - this is 
because the resistivity correction is linear in <j£, , which 
is modified for small-angle scattering [see Eq. Q54p] . 

Thus our results for the classical contribution to the 
magnetoresistance can be summarized as follows. For 
weak damping, lo c t* » 1, and 1 <C (qR c ) 2 *C — the 
resistivity correction is: 



S c Pxx 2v 2 f oo 2 tt* c2 ^ J 2 N {qR c ) 



Po 



£2 S T 

7V7 -L 



2 ~*2 



N=- 



+ (lo - Nlo c ) t 



(58) 



For strong damping, lo c t* <C 1, and 1 <C (qR c ) 2 "C — 
S c p x 



Po 



2ql£ z 



1 + 2e cos — sm(2qR c ) 



(59) 
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whilst for lj c t* <§; 1 and (qR c ) 2 S> 



5 c Pxx 
Po 



= 2ql£ 2 



1 + 2e ^ct s C os I — I sin(2(7_R c 

LO r 



(60) 



Summation over ±w and ±q for the imaginary part 
of K ensure that S c p xy and S c p yy vanish as for short 
range scattering [see Eq. (|36fl ]. The three resistivity 
regimes identified above correspond respectively to high 
[Eq. iJSEl], intermediate [Eq. <|53|l]. and low [Eq. ({67$ ] 
magnetic fields. The crossover between high and interme- 
diate magnetic fields is at lo c t* ~ 1, which is equivalent 
to oj c t ~ {qVjz , and the crossover between intermedi- 
ate and low magnetic fields is when (qR c ) 2 ~ — , which 
is equivalent to w c t ~ ql^J^. The number of damped 
oscillations at low magnetic field has recently been cal- 
culated elsewhere^ As for short range scattering, these 
results are for the regime where channeled orbits can be 
ignored, which requires that the Lorentz force should be 
stronger than the screened SAW field, i.e. £ < — — . 

Both Eqs. and reduce to the results found by 
Mirlin and WolflaSI in the limit that u> — > 0, provided one 
replaces lu/s by q in Eq. (|5*51) . To extrapolate this result 
to lower magnetic fields, the best we can do is a simi- 
lar procedure to Mirlin and Wolfle, which is to replace 
e~"c-r' by the damping factor D discussed in Sec. IIIIBl 
leading to Eq. I|60|) . In the limit that B — > this reduces 
to the same behaviour as in the case of isotropic scat- 
tering, and the magnetoresistance correction is given by 
Eq. (|40|l . We neglect the possibility of channeled orbits, 
as discussed in Sec. IIII A 21 
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FIG. 4: Magnetoresistance correction for small-angle scatter- 
ing, with lot — 7, vf/s = 63, and £ = 0.02. The curve is 
determined by splicing the expressions for the resistance in 
the strong and weak damping limits (lo c t* < 1 and lo c t* > 1 
respectively), which are indicated in the figure. 

In Figs. 01 [S] andElwe show the behaviour of the mag- 
netoresistance as determined by splicing the results of 
Eq. H58|) and Eq. I|59|l for the weak and strong damping 



cases respectively. In Fig 0] we use the same sample pa- 
rameters as in Fig. |2 whilst both Figs. and HJ1 are for 
v F /s = 84, with ujt = 60 in Fig.E3and lot = 600 in Fig.© 
for the sample parameters in Ref. these correspond to 
uj = 2n x 10 GHz and 2tt x 100 GHz respectively. Note 
that, owing to our requirement qR c lu c t, o ur results 
should not be trusted in the region lu c /lu < 'v f / '(suit) 
(which is of order 1 for the parameters used), i.e. at low 
magnetic fields. However, our expectation is that mag- 
netoresistance oscillations should be strongly damped in 
this regime, as described in Eq. H60[). We also note that 
for io c I to 3> 1, the resistance correction in the case of 
small- angle scattering is enhanced over that expected for 
isotropic scattering with the same value of the transport 
time (i.e. the same 2DEG mobility). 




FIG. 5: Magnetoresistance correction for small-angle scatter- 
ing, with lot = 60, vf/s — 84, and £ — 0.003. The curve 
is determined by splicing the expressions for the resistance in 
the weak and strong damping limits, which are indicated in 
the figure. 
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FIG. 6: Magnetoresistance correction for small-angle scatter- 
ing, with lot = 600, v F /s = 84, and £ = 0.002. The curve 
is a splice of the formulae for the strong and weak damping 
limits, and the range of validity for each region is indicated. 
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IV. QUANTUM KINETIC CONTRIBUTION TO 
SAW-INDUCED MAGNETORESISTANCE 

The absorption of SAWs by electrons in two dimen- 
sions changes their steady state distribution over energy. 
Characteristics that are independent of energy cause no 
additional changes to magnetoresistance beyond those 
previously described. However, when uj c t 3> 1, Lan- 
dau level quantization is present, leading to an oscillatory 
energy-dependence of the electron DOS for overlapping 
Landau levels of 7(e) = [1 — L cos (2ire/hui c )) 7, where 
r = 2e '* c ' r « <C 1, and T q is the quantum lifetime of the 
Landau levels42 [A calculation of the density of states 
for a long-range disorder potentials^ shows that T q can 
be field-dependent.] This gives rise to a quantum contri- 
bution to the geometrical commensurability oscillations 
which can persist up to high temperatures. The oscilla- 
tions in the DOS impose oscillations on the electron elas- 
tic scattering rate, T _1 (e) = r _1 7(e)/7, which in turn 
gives a contribution to the observable conductivity. At 
low temperatures, fcgT < hcu c , the DOS oscillations lead 
to Shubnikov-de Haas (SdH) oscillations in conductivity. 
At high temperatures fcflT > hio c , thermal broadening 
smears out the SdH oscillations, but the quantum con- 
tribution can remain as a non-linear effect after energy 
averaging. 

We follow a similar approach to Dmitriev et 
to study the quantum kinetic equation to obtain the 
first quantum correction to the classical magnetoresis- 
tance induced by SAWs. This change in the distribu- 
tion is oscillatory in energy and leads to a contribution 
to the DC conductivity that oscillates as a function of 
u)/u> c . The effect depends on the efficiency of energy 
relaxation and hence is proportional to Ti„; this effect 
dominates the quantum corrections discussed in Refs.lT^I 
when Tin/rq 3> 1. If this is the case, the energy relax- 
ation time is long compared to the Landau level lifetime, 
allowing a strongly non-equilibrium electron distribution 
as a function of energy to arise. 

The photoconductivity cr p h determines the longitudi- 
nal current flowing in response to a DC electric field in 
the presence of SAWs: 



where 



13.14 



En 



CFt, 



ph|Eoo| , 



(61) 



and can be related to the photoresistivity by p p h — 
PxyCpht where p xy = eB/n e . To calculate the photocon- 
ductivity, one integrates over the distribution function 



o-ph = 2 J dea 00 (e) [-d e f(e)} 
and in the leading approximation, 



(62) 



^oo(e) = cr 00 — — , 



(63) 



7 00 



2 2 

e 7 v f 



(64) 



is the Drude conductivity. In obtaining a solution to 
the problem, we are only interested in effects due to the 
non-trivial energy dependence of the electron distribution 
function /(e). 

To perform this calculation we need the classical dy- 
namical conductivity, from which we consider the energy 
dependence of the density of states 7 and the momentum 
relaxation time r. We calculate the classical dynamical 
conductivity below in Sec. lIV A1 and then use it to derive 
the quantum energy balance equation and magnetoresis- 
tance correction in Sec. IIVBI 



A. Classical dynamical conductivity 

The magnetic field dependence of the resistivity change 
reflects the form of the SAW attenuation by the 2D elec- 
trons determined by the real part of the longitudinal dy- 
namical conductivity <j wg , 



Re {(Juq} = is 2 re 2 Re 



K 



1 - K 



(65) 



We obtain a uq from considering attenuation in the form 



@Ljq\Eijq\ , 



and we use the expression in Eq. (|65|) in the following 
subsection. In analogy with our previous work, we find 
that the equation for a^q for small-angle scattering is 

Re{cr wg } = 7s 2 r*e 2 Re|A v | . (66) 



B. Quantum energy balance equation 

The quantum energy balance equation in Ref. ^| is 
stated without detailed derivation. Here we provide a 
simple derivation as an alternative to the approach used 
in Ref. [Lj If we consider the energy balance in a 2DEG 
due to absorption and emission of SAWs, then the rate 
at which energy is added is Q(t) = j uq (t) ■ E uq (t), and 
3uq(t) = o"wqE wg (i), which when we average over time 
gives Q = 5 0^,5 1 E wg 1 2 We can express the rate of absorp- 
tion or emission processes involving energy levels i and j 
as 



2tt 



(67) 



using Fermi's golden rule, where My is the matrix ele- 
ment between states i and j. We assume that My = Al 
and calculate the total power absorption for the classical 
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case (7 constant), which we equate to Q = ^<7^ q \E uq \ 
Hence 



2tt 



-\M\ 2 tkJj(tkJj), 



(68) 



where one factor of ftio is the energy added per photon 
absorbed, 7 is the final density of states, and hwy is the 
number of electrons below the Fermi level that are avail- 
able to make a transition. This gives a formula for \M\ 2 : 



' M ' 2 " si*?- (69 » 

We now consider the change in the number of electrons 
with energy e — -> e + de allowing the density of states to 
be energy dependent, 7(e): 

= Y \ M \ 2 IE /( £ ± ± M*7(e)(l - f(e)) 

-/(e)7(e)de7(e ± ^)(1 - /(e ± M) j 
-^)-Me)- (e)fc 

(70) 

When we require that the distribution be stationary, the 
left hand side of the equation vanishes and when we di- 
vide by 7(e) we get 



* w,|2 *Hr E 7(«±M ± M - « - I/(£) " /r(£)] - 



2S 2 w 2 7 2 

(71) 

If we want to include the effects of a DC field as well, 
then we can see the form of the DC term from the oj — > 
limit of the right hand side of Eq. (|70|l , and we get as a 
final result the balance equation from Ref. H3 . 



2ft 2 w 2 7 2 

|ffoo|Voo d 
7(e) 9e 



]T 7 (e±M [/(e±M-/(e)] 



7(e) 2 d 
7 2 <9e 



[/(e) -Me)] 



(72) 



We could also have obtained the same equation by treat- 
ing a uq as energy dependent and then pulling out fac- 
tors of the energy dependent density of states from 
7 (e) and r _1 (e + hoj), which is equivalent to the ap- 
proach in Ref. Ila . This second approach works provided 
\N — lo/lu c \ < 1/lo c t (i.e. we can ignore the 1 in the de- 
nominator of the sum over N harmonics in cr^ q ). At fre- 
quencies closer to resonance than this, we can no longer 
ignore the 1 in the denominator and the classical ex- 
pression lacks accuracy. However, for large uj c t this is a 



relatively small region of magnetic fields, and is not rel- 
evant at the magnetic fields at which ZRS mimima are 
observed. 

We use the expression for a uq from Eq. I|t)5l) (short- 
range disorder potential) or Eq. I|66l) (long-range dis- 
order potential), which contains the geometric and fre- 
quency commensurability oscillations, and <Tq we know 
as the Drude resistivity [Eq. Introduce the follow- 

ing quantities 



V 



Q = n 



2a^ q \E uq \ 2 
47rV D |£ o| 2 



?i 2 w 2 7 

We are then left with the following equation to solve 



(73) 
(74) 



V 

4 



El 1 - 11 

± v 
47r 2 7 (e) de 



<^^^H(/(e±M-/(e)) 



ftw r 



/27re\\ 2 d , . 



/(e) -/T(e). 



(75) 



Now, let/(e) = / (e)+r/i(e)cos 
M<0, *> /(e) c f T (e) + f osc (e) 



where / (e) ~ 
and 



hiu c 

When we do this, 
then sum over v and neglect all derivatives of /1 that arise 
from Taylor expanding f(e ± Huj) (since this is assumed 
to be a smooth function on the energy scale of ksT), 
and retain only lowest order terms in T, we get a linear 
equation in / osc which is trivial to solve and leads to 
(rioting that Offj- is also ignored) an identical result to 
Ref. " 



mg 1 

M 



ft* r Bfr . ( 2« \ ^Bin(^) + 4Q 
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We then use the expression for the photoconductivity, 
Eq. (|r]2*)l and get the isotropic change in the DC conduc- 
tivity 
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(77) 



Note that in the approach outlined above we made 
no detailed assumptions about the momentum relaxation 
apart from n n 3> r, such that the momentum relaxation 
is efficient in making the distribution isotropic. Thus 
the only difference in the expression for a short or long- 
ranged potential is in the form of a ujqi which enters in 
the expression for V in Eq. i|73|) . We now present numer- 
ical calculations of the resistivity change derived above in 
the presence of SAWs for both isotropic and small-angle 
scattering. All these calculations are for the case Q = 0. 
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C. Numerical results 

In this section we show the resistivity as it is modified 
by the quantum correction alone. We discuss the situ- 
ation where both classical and quantum corrections are 
important in Sec. IV Al 



1. Isotropic scattering 

We illustrate numerical results for the magnetoresis- 
tance oscillations due to the combination of SAWs and 
density of states modulations in Fig. We use Eq. f2"%|l 
for K and do not ignore the 1 in the denominator of 
the summand. For large ujt, this should give accurate 
results except for a very small range of magnetic fields 
around each of the resonances where uj/uj c equals an in- 
teger. Note that we divide the expression in Eq. (|77Jl by 
1 + r 2 /2, to compare the resistivity with and without 
SAWs, assuming that there are density of states modula- 
tions in both cases. To connect to our previous notation, 
we note that 

V = 2 (^) 2 ( W r)( W r m )£ 2 Re . (78) 

In Fig. we show the change in the magnetoresis- 
tance for ujT q — 6 and ujr in — 302 (usi ng the value of 
Ti n = 7.64 x 10 -10 s estimated in Ref. |lj) for the param- 
eters that we calculated the classical magnetoresistance 
corrections displayed in Figs. [21 and As for the case of 
a microwave field, there are peaks in the resistance when 
lu/u> c is close to an integer, but with Sp xx — when u>/uj c 
is exactly equal to an integer. The peaks (and dips) are 
modulated by geometric commensurability oscillations, 
and this leads to a new class of ZRS, of the type pre- 
dicted for to <C lo c in Ref.l28l For lot = 7, vp/s = 63, the 
ZRS occur for £ > 0.0015 and for lot = 60, v F /s = 84, 
ZRS occur for £ > 0.001. 

At large lo c /lo ^> 1 (not shown) the resistance change is 
negative as we predicted analytically previously. 28 How- 
ever, the magnitude of the resistance change is very much 
smaller than when uj/lj c > 1, so a very large value of £ 
would be required to attain ZRS at uj c /lj ^ 1. 



2. Small angle scattering 

In Fig. |S| we show data for the quantum magnetoresis- 
tance correction when there is small-angle scattering for 
the same sample parameters and frequencies as in Figs. 01 
and[Sl Similarly to the short-range case, we can relate V 
to our previous notation: 

V = 2 (g) 2 (LOT*)^T tn )£ 2 Re{k} , (79) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

co c /co 

FIG. 7: Quantum contribution to the magnetoresistance for 
ujt = 7, vf/s — 63 and cut — 60, vf/s = 84, both with 
uiT q = 6, £ = 0.0015, and uTi n = 302. 



where the change from Eq. I|78|l reflects that we use a 
different cr wg in the two cases. 
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FIG. 8: Quantum contribution to magnetoresistance for ujt = 
7, vf/s = 63 and for ujt — 60, vf/s = 84. In both cases 
uiT q = 6, £ = 0.003, and ujT in = 302. The arrows indicate the 
range of validity of the results. 



If lo c /lo < 1 then there are ZRS of the type that are 
usually observed with microwaves, however there are an- 
other class of magnetoresistance oscillations at larger val- 
ues of Ld c /uj, which are of the sort we predicted in Ref.l28l 
and are much less strongly damped than in the isotropic 
scattering case. We note that for ujt = 60, vp/s = 84, 
the threshold value for ZRS at uj c /lu > 1 is £ ~ 0.003. 
In Fig. |Slwe clearly demarcate the values of lu c /lu where 
our analytic formulae are valid {qR c /u) c T 1). 
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V. DISCUSSION 
A. Combination of classical and quantum effects 

At low frequencies (oj <C oj c ) we can use the approach 
in Sec. IIVBI to obtain an analytic expression for the 
quantum resistivity correction, starting with the low fre- 
quency limit of the balance equation for vanishing dc 
field. This calculation was already considered in Ref. |2^| 
and we quote the result here for completeness. We find 
a non-vanishing addition to both diagonal components 
of the conductivity (and, therefore, also of the resistiv- 
ity) even when ksT 3> hoj c . This generates isotropic 
magneto-oscillations for short-range disorder 



S q p 
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5 q <J aa 



T 



AitTe F 



huj r 



£ 2 J 2 (qR c ), (80) 



(where Co = 1/po) m addition to the anisotropic classical 
commensur ability effect. Together, they yield the result: 
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(81) 



valid when r _1 < uj <C u> c , where v = 2ef/Soj c is the 
filling factor. For a long-range potential, the oscillations 
are very similar, for r* _1 < oj <C w CJ with t* replacing r 



to give 
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The actual magnetoresistance trace observed in an exper- 
iment will have both anisotropic classical and isotropic 
quantum contributions. The parameter that controls the 
overall behaviour for both isotropic scattering and small- 
angle scattering is 



7} = (2ttTvs/v f ) y/nn/r. 



(83) 



We can use r\ to identify three different magnetoresistance 
regimes. If 77 < 1, the observed change in the resistivity 
is in the direction of SAW propagation. If r\ ~ 1, then the 
observed change in the resistivity is in the direction per- 
pendicular to the SAW wavevector, and if 77 3> 1, then 
the resistance correction is isotropic and negative. For 
the parameter values used in Figs. and |SJ the param- 
eter rj is large compared to one, such that the quantum 
contribution dominates the classical corrections. 

The quantum contribution to the resistivity that we 
discuss in Sec. II VI is not the only possible contribution 
to the resistivity that can arise from quantum effects ne- 
glected in a classical calculation. A static periodic poten- 
tial affects the equilibrium density of states, which can 



lead to corrections to p yy as well as to p X x 20]42 These 
corrections were ignored in our calculations, but at least 
in the static case appear to be at most the order of mag- 
nitude of the classical contribution. Hence the quantum 
effect we discuss here should dominate that quantum cor- 
rection for 77 > 1, provided r» n 3> t and v > vf/2tts. 



B. Experimental Implications 

There are four situations we have discussed in this pa- 
per: the combinations of either short- or long-range dis- 
order and quantum and classical corrections to the mag- 
netoresistance due to SAWs. We first focus on the case 
of isotropic scattering. The most pronounced geometric 
and temporal oscillations in the magnetoresistance are 
in the region lo c /uj < 1 regardless of whether the oscil- 
lations are quantum or classical in origin. The classi- 
cal effects dominate for short inelastic scattering times 
and lead to anisotropic, positive magnetoresistance cor- 
rections, whilst quantum effects dominate for large in- 
elastic scattering times, and can lead to ZRS which are 
modulated by geometric commensurability oscillations. 
For small-angle scattering, geometric commensurability 
oscillations in the regime uj c /oj < 1 are very strongly 
damped in both the classical and quantum cases, unless 
(qR c ) 2 /lu c t < 1 in this region. This condition is not met 
in current high quality 2DEGs in which small-angle scat- 
tering dominates. However, it is currently possible to 
achieve (qR c ) 2 / uj c t < 1 for lo c /uo > 1, and these geomet- 
ric commensurability oscillations should be observable ei- 
ther in the quantum or classical regimes when lo c /lo > 1. 
Similarly to isotropic scattering, quantum effects domi- 
nate for large inelastic scattering times (compared to the 
elastic scattering time) and classical effects are more im- 
portant for short inelastic scattering times. It appears 
that the ZRS at low frequencies we predicted in Ref. 
for isotropic scattering, are hard to achieve if the scatter- 
ing is isotropic, but should be much more readily achiev- 
able in samples where small-angle scattering dominates, 
which is the experimentally relevant situation. 

One theoretical prediction for ZRS is that they re- 
quire inhomogeneous current flow in a 2DEG,-- 1 which 
appears to have recently been observed for microwaves^ 
A large enough SAW- induced change \Sa xx \ > &o result- 
ing in negative local conductivity would also require the 
formation of electric field/Hall current domains. Since 
the anisotropy in Eq. ijST} suggests that such conditions 
can be achieved most easily in the conductivity compo- 
nent along the SAW wavevector, we expect that domains 
would form with current flowing perpendicular to the di- 
rection of SAW propagation, and their stability would 
depend on the sample geometry. For a SAW with the 
wavevector directed across the axis of a Hall bar, current 
domains can be stabilized by ending in ohmic contacts. 
For a wave propagating (or standing) along the Hall bar, 
current domains would have to orient across the bar di- 
rection and terminate at the sample edges (destabiliz- 
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ing them), leading to a finite resistance. Finally, the 
anisotropy would not support a zero-conductance regime 
in a Corbino geometry. 

In addition to the parameter ranges that are optimal 
for observing SAW- induced ZRS, there are several other 
experimental issues we would like to mention. Firstly 
there is the observation by Kukushkin et al& that mi- 
crowave irradiation of high quality 2DEGs can lead to 
SAW generation. Thus one might want to consider the 
effects of microwaves and SAW simultaneously - such a 
scheme might also allow for probing a 2DEG at frequen- 
cies other than qs for a given SAW wavenumber. How- 
ever, this calculation is beyond the scope of the present 
work4£ 

If increases in SAW frequency and 2DEG mobility al- 
low lot 3> (qRc) 2 for lo/lo c > 1, the geometric oscillatory 
structure similar to that observed for isotropic scattering 
i.e. peaks when lv/lv c is an integer should be visible in 
both classical and quantum contributions to magnetore- 
sistance. This is currently not yet observable for samples 
in which small-angle scattering dominates (as appears 
to be the case for the samples used in Refs. 00 where 
ZRS were observed) because lu c t* 3> 1 when lo/lo c > 1. 
Either samples in which isotropic scattering dominates 
are required to see these effects, or higher quality sam- 
ples in which (qR c ) 2 /uj c t 1 when lo/lo c > 1 are re- 
quired. In terms of trying to improve the possibility of 
observing interesting SAW induced magnetoresistance ef- 
fects in the frequency range lo > lo c in samples where 

small-angle scattering dominates, the following consider- 

i 

ations may be helpful. Since vp/s oc n| and lot oc lo^i 
(and it is found that [i oc n® m the highest quality 
GaAs/AlGaAs 2DEGsii), the best hope for observing ef- 
fects with lot 3> vp/s appears to be by achieving higher 
SAW frequencies. 

C. Conclusions 

In conclusion, we have demonstrated a new class of 
magnetoresistance oscillations caused in a 2DEG by 



SAWs. We have shown that at w < w c the effect con- 
sists of contributions with competing signs: (i) a classi- 
cal geometric commensurability effect analogous to that 
found in static systems with positive sign, and (ii) a quan- 
tum correction, with negative sign for either isotropic or 
small-angle scattering. The latter result suggests that 
SAW propagation through a high mobility electron gas 
may generate a sequence of zero-resistance states (ZRS) 
linked to the maxima of Jq (qR c ) for strong enough SAW 
fields. We find that in the region lo <C lo c , SAW-induced 
ZRS states are much more likely to be observed in sam- 
ples for which small-angle scattering dominates isotropic 
scattering. Whilst this prediction concerns the low- 
frequency domain lo < lo c , such ZRS would be formed 
via the same mechanismAi*ii as the microwave-induced 
ZRS at lo > lo c . In the regime lo > lo Ci we find that there 
are geometric oscillations superposed on ZRS if there is 
isotropic scattering. If there is small-angle scattering, 
such oscillations are unlikely to be seen in 2DEGs at the 
present time. Hence the optimal parameter region to 
search for SAW induced ZRS in 2DEGs with long-range 
disorder which show geometric modulation is for lo < lo c 
and (qR c ) 2 /lo c t < 1. 
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